Campbell diagrams of weakly anisotropic 

flexible rotors 



By O.N. KiRiLLOvf 

We consider an axi-symmetric flexible rotor perturbed by dissipative, conser- 
vative, and non-conservative positional forces originated at the contact with the 
anisotropic stator. The Campbell diagram of the unperturbed system is a mesh-like 
structure in the frequency-speed plane with double eigenfrequencies at the nodes. 
The diagram is convenient for the analysis of the traveling waves in the rotating 
elastic continuum. Computing sensitivities of the doublets we find that at every par- 
ticular node the unfolding of the mesh into the branches of complex eigenvalues in 
the first approximation is generically determined by only four 2x2 sub-blocks of the 
perturbing matrix. Selection of the unstable modes that cause self-excited vibrations 
in the subcritical speed range, is governed by the exceptional points at the corners 
of the singular eigenvalue surfaces — 'double coffee filter' and 'viaduct' — which are 
sharply associated with the crossings of the unperturbed Campbell diagram with 
the definite symplectic (Krein) signature. The singularities connect the problems of 
wave propagation in the rotating continua with that of electromagnetic and acoustic 
wave propagation in non-rotating anisotropic chiral media. As mechanical examples 
a model of a rotating shaft with two degrees of freedom and a continuous model of 
a rotating circular string passing through the eyelet are studied in detail. 

Keywords: Campbell diagram; flexible rotor; dissipation-induced instabilities; 
subcritical flutter; symplectic (Krein) signature; non-Hermitian degeneracies 



1. Introduction 

Bending waves propagate in the circumferential direction of an elastic body of rev- 
olution rotating about its axis of symmetry (Bryan 1890, Southwell 1921, Shrini- 
vasan and Lauterbach 1971, Genta 2007). The frequencies of the waves plotted 
against the rotational speed are referred to as the Campbell diagram (Campbell 
1924, Genta 2007). The spectrum of a perfect rotationally symmetric rotor at 
standstill has infinitely many double semi-simple eigenvalues — the doublet modes. 

Indeed, for M^^^ 9 A = diag(cj?, cjn) and R = f sinO \ ^^^^ restriction 

^ ^' \ -sm0 cos6' / 

R^AR = A imposed by equivariance of the equations of motion with respect to 

the action of the circle group implies = see, e.g., (Dellnitz et al. 1992, Na- 

gata and Namachchivaya 1998). By this reason the Campbell diagram contains the 

eigenvalue branches originated after splitting of the doublets by gyroscopic forces 

(Bryan 1890). The branches correspond to simple pure imaginary eigenvalues and 

intersect each other forming a spectral mesh (Giinther and Kirillov 2006) in the 
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Figure 1. (a) The Campbell diagram of the unperturbed system (2.2) with 6 d.o.f. in case 
oi LUi = 1, UJ2 = and uo^ = 6; (b) the Campbell diagram and (c) decay rate plots for 
the stiffness modification kKi with k = 0.2; (d) unfolding the Campbell diagram due to 
perturbation with the matrices K = Ki, D = Di, N = Ni and k = 0.2, 5 = O.f, and 
V — 0.2, (e) the corresponding Campbell diagram and (f) decay rate plots. 



frequency-speed plane with the doublets at the nodes. Fig. 1(a). Perturbations of 
the axially symmetric rotor by dissipative, conservative, and non-conservative posi- 
tional forces, caused by its contact with the anisotropic stator, generically unfold the 
spectral mesh of pure imaginary eigenvalues of the Campbell diagram into separate 
branches of complex eigenvalues in the (il, ImA, ReA)-space, see Fig. 1(d). Never- 
theless, the eigenvalue branches in the perturbed Campbell diagram can both avoid 
crossings and cross each other. Fig. 1(e). Moreover, the real parts of the perturbed 
eigenvalues plotted against the rotational speed — decay rate plots (Genta 2007) — 
can also intersect each other and inflate into 'bubbles'. Fig. 1(f). This complicated 
behavior is difficult to predict and even to interpret according to the studies of nu- 
merous mechanical systems (Ono et al. 1991, Chen and Bogy 1992, Mottershead and 
Chan 1995, Yang and Hutton 1995, Tian and Hutton 1999, Xiong et al. 2002, Canchi 
and Parker 2006, Young and Lin 2006, Genta 2007, Lesaffre et al. 2007, Spelsberg- 
Korspeter et al. 2009). The present work reveals that the unfolding of the Campbell 
diagrams is determined by a limited number of local scenarios for eigenvalues as a 
function of parameters, which form stratified manifolds. 



2. A model of a weakly anisotropic rotor system 

In general, the imperfections in the rotor and stator complicate the linearized equa- 
tions of motion making them non-self-adjoint with time-dependent coefficients. Nev- 
ertheless, an axially symmetric rotor with an anisotropic stator as well as an asym- 
metric rotor with an isotropic stator are autonomous non-conservative gyroscopic 
systems (Genta 2007). Neglecting the centrifugal stiffness without loss of generality. 
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we consider the finite-dimensional anisotropic rotor system 

X + {2nG + (5D)x + (P + r^^G^ + kK + z/N)x = 0, (2.1) 

which is a perturbation of the isotropic one (Nagata and Namachchivaya 1998) 

x-|-2nGx+(P-|-Q^G^)x = 0, (2.2) 

where x = K^", P=diag(a;i , wj, w|, . . . , w^, w^) is the stiffness matrix, and G = 
— G^ is the matrix of gyroscopic forces defined as 

G = blockdiag(J, 23,..., nJ), J = ( J ~J ) • (2-3) 

The matrices of non-Hamiltonian perturbation corresponding to velocity-dependent 

dissipativc forces, D = D^, and non-conscrvativc positional forces, N = — N"^, 
as well as the matrix K = of the Hamiltonian perturbation that breaks the 
rotational symmetry, can depend on the rotational speed fl. The intensity of the 
perturbation is controlled by the parameters S, n, and v. Putting k = and v = Q 
in (2.1) yields the model considered in (Nagata and Namachchivaya 1998). 

At O = the eigenvalues iiwg, Wg > 0, of the isotropic rotor (2.2) are double 
scmi-simplc with two linearly independent eigenvectors. The sequence of the fre- 
quencies ujs, where s is an integer index, is usually different for various bodies of 
revolution. For example, ojg = s corresponds to the natural frequency /« = 

of a circular string of radius r, circumferential tension P, and mass density p per 
unit length (Yang and Hutton 1995). 

Substituting x = uexp(At) into (2.2), we arrive at the eigenvalue problem 

Lo(^])u := (lA^ -h 217GA -h P -h Vl^Q,'^)\x = 0. (2.4) 

The eigenvalues of the operator Lq are found in the explicit form 

Xf = iu)s + isQ, X7 = —ius + isfl, X~ = iiOs — isfl, Xf = —itOs — is^, (2.5) 

where the overbar denotes complex conjugate. The eigenvectors of Aj" and Xs are 

u+ = (-j, 1,0,0,..., 0,0)^, u+ = (0,0,...,0,0,-i,l)^ (2.6) 

where the imaginary unit holds the (2s — l)st position in the vector u+. The eigen- 
vectors, corresponding to the eigenvalues A J and Xf , are simply uj = uf 

For > 0, simple eigenvalues A+ and A^" correspond to the forward and back- 
ward traveling waves, respectively, that propagate in the circumferential direction 
of the rotor. At the angular velocity Qf~ = ujs/s the frequency of the ,sth backward 
traveling wave vanishes to zero, so that the wave remains stationary in the non- 
rotating frame. We assume further in the text that the sequence of the doublets 
ius has the property tUs+i — > ^ which implies the existence of the minimal 
critical speed Qcr = — '^i- When the speed of rotation exceeds the critical 
speed, some backward waves, corresponding to the eigenvalues XJ , travel slower 
than the disc rotation speed and appear to be traveling forward (reflected waves). 
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In Fig. 1(a) the spectral mesh (2.5) is shown for the 6 d.o.f.-system (2.2) with 
the frequencies uj\ = 1, uj2 = and = 6 that imitate the distribution of the 
doublets of a circular ring (Canchi and Parker 2006). To illustrate typical unfolding 
of the Campbell diagram, we plot in Fig. l(d)-(f) the eigenvalues of the 6 d.o.f.- 
system (2.1) with K = 0.2, 5 = 0.1, v = 0.2, wi = 1, ^2 = 3, and = 6 for the 
specific symmetry-breaking matrix K = Ki, (Ki = K^), whose non-zero entries 
are kn = 1, ki2 = 2, fcis = 1, ku = 2, ^22 = 1, fes = 3, k2A = 4, /C33 = -3, 
kii = —2.5, = 4, = 2, and for the matrices D = Di and N = Ni, where 

1 1 -1 -3 8 \ 
2 3 2 4 

2 7 1 3 

3- 7 8 2 
2-1-8 2 

4- 3-2-20/ 

In the following wc classify and interpret the typical behavior of the eigenvalues 
of weakly anisotropic rotor system (2.1) with the use of the perturbation formula 
for the doublets of the spectral mesh (2.5), which we derive in the next section. 
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3. Perturbation of the doublets 



Introducing the indices a,/?, e, cr ~ ±1 we find that the eigenvalue branches = 
iaijJs + ies^ and = ifiujt + iatQ cross each other at f2 = fio with the origination 
of the double eigenvalue Aq = iu)o with two linearly-independent eigenvectors 



and u^, where 



Wo 



aaust — (3eujts 



(3.1) 



at — es at — ss 

Let M be one of the matrices D, K, or N. In the following, we decompose the 
matrix M e M2nx2n j^^^ ^2 yocks M^t € M2x2^ wheve s,t = 1,2, ... ,n 



M = 



( * 



M,, 



* 



"^2s-l,2t-l ?Tl2s-l,2t 
m2s,2t-l m2s,2t 



(3.2) 



Note that D,t = Df,, K,* = Kf,, and N«t = -Nf,. 

We consider a general perturbation of the matrix operator of the isotropic rotor 
Lo(fi) -I- AL(ri). The size of the perturbation AL(0) = JAD -|- kK -|- i/N ~ e is 
small, where e = ||AL(rio)|| is the Frobenius norm of the perturbation at Q = Oq- 
For small Afi = jfi — fiol and e the increment to the doublet Aq = ioj^) with the 
eigenvectors w% and w", is given by the formula det(R-|- (A — Ao)Q) = (Kirillov 
et al. 2005, Kirillov 2008), where the entries of the 2x2 matrices Q and R are 



Q 



i£(T 
St 



T>ea 



= (2ia;o(u' 



^<+20o(uD^ 
)'^G<+2f2o(u; 



(3.3) 
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Calculating the coefficients (3.3) with the eigenvectors (2.6) we find the real and 
imaginary parts of the sensitivity of the doublet Aq = icoo at the crossing (3.1) 



ReA = _if^^^!^ + ^^^±./M-Rec 
8 \ atOs pujt 



T ^ ^^^t , ^ , , ^ tiK,, , tYKtt\ ^ / |c|+Rec 

ImA = -o + ^(as + ia) + -^^ + -^J±y^— , (3.4) 

where c = Rec + ilmc with 

Imc = [se — t(T)A\l 

SuJsiOt 

(gg^gtrKit - Po JttrKss){a(^siToaBi - PoJtlmAi) 



2 



, ta- se PtOstrKtt - awttrK, 



+ Ctp — K 



(awJmSi - putlmAif + Aaf3usUJt{{ReA2y + (ReSs)^) 



64a;2o;2 



(3.5) 



The coefficients Ai , A2 and Bi , B2 depend only on those entries of the matrices 
D, K, and N that belong to the four 2x2 blocks (3.2) with the indices s and t 

Ai = (5AotrDss+KtrKj5s+e2?:;/n2s-i,2s, 

A2 = <ji/trNsfIe<^+z(i5AotrDstJ£cr+KtrKsf Jscr), 

Bi = SXotT'Dtt+KtiKtt+cT2ii'n2t-i,2t, 

B2 = autr'NstJscT-i{S>^otr'Dsd£a+KtrKstIea), (3.6) 

where 



Therefore, we have identified the elements of the perturbing matrices that control 
practically important eigenvalue assignment (Ouyang 2008) near every particular 
node {flo,u}o) of the spectral mesh. 



4. MacKay's eigenvalue cones and instability bubbles 

Modification of the stiffness matrix induced by the elastic support interacting with 
the rotating continua is typical in the models of rotating shafts (Shieh and Masur 
1968); computer disc drives (Ono ct al. 1991, Chen and Bogy 1992), circular saws 
(Yang and Hutton 1995, Tian and Hutton 1999, Xiong et al. 2002), car brakes 
(Mottershead and Chan 1995, Herve et al 2008, Spelsberg-Korspeter et al. 2009), 
and turbine wheels (Genta 2007, Lesaffre et al. 2007). 
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Figure 2. Eigenvalue surfaces (MacKay, 1986) and (bold lines) their cross-sections in the 
plane k = const (grey): (a) a near- vertically oriented cone ImA(f2, k) in the subcritical 
range (ReA = 0); (b) imaginary parts forming a near-horizontally oriented cone (4.1) with 
the attached membrane (4.2) and (c) the real parts forming a near-horizontally oriented 



cone (ReA) 



-Rec with the attached membrane ReA = in the supercritical range. 



Assuming (5 = and v = {) in (3.4) wc find that the eigenvalues of the system 
(2.5) with the stiffness modification kK either are pure imaginary (ReA — 0) and 
form a conical surface in the (51, k, lmA)-space with the apex at the point (fio, 0, wq) 



1mA — 0^0 — — I h 



8 V aujs 



n-n 



{se + ta) I = Rec 



(4.1) 



see Fig. 2(a), or they are complex and in the (il, k, ReA)-space their real parts 
originate a cone (ReA)^ = —Rec with the apex at the point (iloiO,0), Fig. 2(c). In 
the (Jl, K, ImA)-space the corresponding imaginary parts belong to the plane 



ImA = ljq 



K ( trKo 



trK» 



(se + ta), 



(4.2) 



which is attached to the cone (4.1) as shown in Fig. 2(b). 

The existence of eigenvalues with ReA depends on the sign of a(3. It is 
negative only if the crossing in the Campbell diagram is formed by the eigenvalue 
branch of the reflected wave and by that of either forward- or backward traveling 
wave. Otherwise, a(3 > 0. Due to the property lOs+i ~ lOs > flf the crossings of 
the reflected wave with the forward- and backward traveling waves occur only in 
the supercritical speed range |r2| > ilcr- The crossings with a/3 > are situated in 
both the super- and subcritical < Qcr) ranges. Therefore, the eigenvalues with 
ReA =/= originate only near the supercritical crossings of the eigenvalue branches 
XI and with a/3 < 0, when the parameters in the [Q, K)-plane are in the sector 
Rec < bounded by the straight lines Rec = 

A{ss ~ ta){n - no) 



i+k2 



s-l+k2 



'-±2 



I (£fc23-l,2t-l+a'fc23.2t)^ + (£fc2s-l,2f-g-fc2»,2t-l)^ 

(4.3) 

Since for a/3 < the cones of the real parts (ReA)^ — —Rec are near-horizontally 
oriented and extended along the K-axis in the (fi, k, ReA)-space, their cross-sections 
by the planes k = const are ellipses, as shown in Fig. 1(c) and in Fig. 2(c). Since a 
part of the ellipse corresponds to the eigenvalues with positive real parts, the ellipse 



Article submitted to Royal Society 



Campbell diagrams of weakly anisotropic flexible rotors 7 

is called the bubble of instability (MacKay and Saffman 1986). Equation (4.3) is, 
therefore, a linear approximation to the boundary of the domain of instability, 
which is divergence (parametric resonance) for Qq = Clf and flutter (combination 
resonance) otherwise. The near-horizontal orientation of the corresponding cones of 
imaginary parts (4.1) in the {fl, k, ImA)-space explains deformation in the presence 
of the perturbation kK of the crossings with af3 < into the branches of a hyperbola 
connected by a straight line in the Campbell diagram, see Fig. 1(b) and Fig. 2(b). 

Near the crossings with a/3 > the perturbed eigenvalues are pure imaginary 
(stability). The corresponding cones of imaginary parts (4.1) are near- vertically 
oriented in the (f2, k, ImA)-space, Fig. 2(a). In the plane k = const this yields the 
avoided crossing (MacKay 1986, MacKay and Saffman 1986), which is approxi- 
mated by a hyperbola shown by the bold lines in Fig. 2(a) (cf. Fig. 1(b)). 

The conical singularities of the eigenvalue surfaces in the Hamiltonian systems 
are traced back to the works of Hamilton himself, who predicted the effect of conical 
refraction of light in birefringent crystals (Hamilton 1833, Berry and Jeffrey 2007). 
Later on, the conical singularities of eigenvalue surfaces were found in atomic, nu- 
clear, and molecular physics (Von Neumann and Wigner 1929, Teller 1937, Mon- 
dragon and Hernandez 1993). Nowadays they bear a name of the Hamilton's diaboli- 
cal points (Berry and Jeffrey 2007). The existence of the two different orientations of 
the eigenvalue cones in the Hamiltonian systems was established in (MacKay 1986). 
This result is based on the works of Williamson (1936) and Krein (1983), who in- 
troduced the signature of eigenvalues known as the symplectic signature in the 
Hamiltonian mechanics (MacKay and Sepulchre 1998) and as the Krein signature 
in a broader context of the theory of Krein spaces (Kirillov et al. 2009). 

To evaluate the symplectic signatures, we reduce (2.2) to y = Ay, where 

A = ( ^^p _0G ) = J2nA^J2„, J2n = ( ) , Y = ( ^ +^Gx ) " 

(4.4) 

The Hamiltonian symmetry of the matrix A impHes its self-adjointness in a Krein 

—T 

space with the indefinite inner product [a, b] = b 32n^, a, b e C^". The matrix 
A has the eigenvalues given by the formulas (2.5) with the eigenvectors 

a++ ^ a+- = ( 

A^u+ + nGut J' ^ { Xjuj + nGuJ 

where the vectors are determined by expressions (2.6). Since i[a++,a++] = 
i[aj ,a+~] = 4:Uig > 0, the eigenvalues A+ and A~ of the forward and backward 
traveling waves acquire positive symplectic (Krein) signature. The eigenvalues A^ 
and A7 of the reflected waves with ?'[a7+,a7+] = i[a7^,a^ ] = —4ujs < 0, have 
the opposite, negative symplectic (Krein) signature (MacKay 1986, MacKay and 
Sepulchre 1998). The signature of an eigenvalue in the Campbell diagram coincides 
with the sign of the doublet at Q = 0, from which it is branched, and does not 
change with the variation of fl. This implies a/3 > and near- vertically oriented 
cones of imaginary parts (4.1) at the crossings of eigenvalue branches with definite 
(positive) signature and af3 < and near-horizontally oriented cones of imaginary 
parts (4.1) at the crossings with mixed signature (MacKay 1986). 

The symplectic signature coincides with the sign of the second derivative of 
the energy, which is a non-degenerate definite quadratic form on the real invariant 
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space associated to a complex conjugate pair of simple pure imaginary non-zero 
eigenvalues (MacKay 1986). Interaction of waves with positive and negative en- 
ergy is a well known mechanism of instability of the moving fluids and plasmas 
(MacKay 1986, Stepanyants and Fabrikant 1989, Hirota and Fukumoto 2008); in 
rotor dynamics this yields flutter in the supercritical speed range, which is known 
as the mass and stiffness instabilities (Mottershead and Chan 1995, Gcnta 2007). 

Therefore, in case when anisotropy of the stator is caused by the stiffness mod- 
ification only, the unfolding of the Campbell diagram is completely described by 
one-parameter slices of the two-parameter MacKay's eigenvalue cones. Since there 
are only two possible spatial orientations of the cones corresponding to either def- 
inite or mixed symplectic signatures, all one has to do to predict the unfolding 
of the Campbell diagram into avoided crossings or into bubbles of instability is 
to calculate the signatures of the appropriate eigenvalues of the isotropic rotor. 
In the following, we develop the MacKay's theory further and show that even in 
the presence of non-Hamiltonian perturbations, all the observed peculiarities of the 
Campbell diagrams and decay rate plots are one-parameter slices of the eigenvalue 
surfaces near a limited number of other singularities whose type is dictated by the 
definiteness of the symplectic signature of the double eigenvalues at the crossings. 



5. Double coffee filter singularity near the crossings with 
definite symplectic (Krein) signature 

Understanding general rules of unfolding the Campbell diagrams of weakly anisotropic 
rotor systems in the presence of dissipative and non-conservative perturbations is 
important for linear stability analysis and for interpretation of numerical data in 
both low- and high-speed applications (Genta 2007). In the latter supercritical flut- 
ter and divergence instabilities are easily excited near the crossings with the mixed 
symplectic signatiirc just by the Hamiltonian perturbations like stiifncss modifi- 
cation. In low-speed applications unfolding of the Campbell diagram is directly 
related to the onset of friction-induced oscillations in brakes, clutches, paper calen- 
ders, and even in musical instruments like the glass harmonica (Spurr 1961, Ouyang 
and Mottershead 2001, Kirillov 2008, Kang et al. 2008, Chevillot et al. 2008, Herve 
et al 2008, Ouyang 2008, Spelsberg-Korspeter et al. 2009). In contrast to the super- 
critical instabilities, excitation of the suhcritical flutter near the crossings with the 
definite symplectic signature by the Hamiltonian perturbations only, is impossible. 
In this case the non-Hamiltonian dissipative and circulatory forces are required for 
destabilization. 

In general, dissipative, ST), and non-conservative positional, vN, perturbations 
unfold the MacKay's eigenvalue cones (4.1) and (ReA)^ = — Rec into the surfaces 
ImA(f7, k) and ReA(f7, k), described by formulas (3.4). The new eigenvalue surfaces 
have singularities at the exceptional points (Keck et al. 2003, Berry and Dennis 
2003). The latter correspond to the double eigenvalues with the Jordan chain that 
born from the parent semi-simple doublet iujo at Q = Qq. In some works numerical 
methods were developed to find the coordinates of these singularities (Jones 1988, 
Skorokhodov 2007). Perturbation of the Hamilton's diabolical points is another 
efiicient way to locate exceptional points (Kirillov et al. 2005). Indeed, condition 
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Figure 3. (a) The 'double coffee filter' singular surface ImA(f2, k) with the exceptional 
points (open circles) and branch cut (bold lines) originated fi'om the MacKay's cone 
(dashed lines) due to mixed dissipative and circulatory perturbation at any crossing with 
the definite symplectic signature; (b) the corresponding 'viaduct' ReA(r2,K). 



; yields their approximate loci in the (fi, K)-plane 



where 



U = 



D 



AuJsUJtita — se) 

ReA2ti-KstJea - ReBztrK.Je 
auJslmBi — PujtimAi 



^EP 



(5.2) 



N 



U^ + af3 



2 

+ a/3 



ReAo 



ReB2 



The crossings with the definite symplectic signature (a/? > 0) always produce a 
pair of the exceptional points. For example, for pure non-conservative {S = 0) and 
pure dissipative {v = 0) perturbation of the doublets at = 0, formulas (5.1) read 



^'■EP,n 



2vn2s 



-1,2s 



''^^•"~^pi(K,,)-p2(K, 



EP,d 



±5 



4s 



^EP.d 



0, 



(5.3) 



where pi_2(Kss) are the eigenvalues of the block K^^ of the matrix K and /ii_2(Dss) 
are those of the block Tigs of D. In case of the mixed symplectic signature (a/3 < 0) 
the two exceptional points exist when N/D > and do not exist otherwise. 

Strong influence of the exceptional points on the stability and their relation 
to the Ziegler's destabilization paradox due to small damping is well recognized 
(Bottema 1956, Kirillov 2007a, Krechetnikov and Marsden 2007, Kirillov 2007b, 
Chevillot et al. 2008). In numerous applications in rotor dynamics (Ono et al. 1991, 
Chen and Bogy 1992, Mottershead and Chan 1995, Yang and Hutton 1995, Tian 
and Hutton 1999, Xiong et al. 2002, Genta 2007, Lesaffre et al. 2007) as weh as 
in hydrodynamics (Or 1991), crystal optics (Berry and Dennis 2003), acoustics 
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(Shuvalov and Scott 2000), and microwave billiards (Keck et al. 2003), the gener- 
alized crossing scenario in the vicinity of the exceptional points has been observed 
(visible also in Fig. l(c,f)) when at the same values of the parameters the imaginary 
parts of the eigenvalues cross, whereas the real parts don't and vice versa. In our 
setting, the conditions for coincidence of imaginary parts of the eigenvalues (3.4) 
arc Imc = and Rcc < and that for coincidence of the real parts arc Imc = 
and Rec > 0. Both real and imaginary parts of the eigenvalues coincide only at 
the two exceptional points (Q^p, fc^p) and (Q^p,fi;^p). The segment of the line 
Imc = connecting the exceptional points is the projection of the branch cut of 
a singular eigenvalue surface ImA(f2, k). The adjacent parts of the line correspond 
to the branch cuts of the singular eigenvalue surface ReA(f2, k). Since simultaneous 
intersection of the different segments of the line Imc = in the (17, K)-plane is not 
possible one observes the generalized crossing scenario (Keck et al. 2003, Kirillov 
et al. 2005) in the planes (Q,ImA) and (f2,ReA) or (k, ImA) and (k, ReA). 

For example, in case of pure non-conservative positional perturbation the real 
parts of the eigenvalues developing near the doublets at fio = cross each other in 
the (O, ReA)-plane at the points of the branch cuts > (k^p„)^ 

ReA = ± '^^sn2s-,, 2s ^ ^ ^^^3^^ ^5 4) 

(Pi(K,,)-P2(K,,))^k2_(^±^^J2 

whereas for < (k^p„)^ they avoid crossing 



At the exceptional points k = K^p „ the (ngc^nvalue branches touch each other 



2 V 

The degenerate crossing (5.6) of the real parts has been observed in the model of a 
rotating circular string (Yang and Hutton 1995, Kirillov 2008). 

Pure dissipative perturbation of the doublets at f^o = yields crossings of the 
real parts at the branch cuts > {^EPd)"^ (ReA, K)-plane and veering of 

the imaginary parts 



ImA = L0,±s^^n^-{n%p^^)'^+O{K), 

ReA = _^!£^± ^ 6k + 0{k^), (5.7) 



^2 



'■EP,d) 

where 7 = 2trKssDss— trKsstrD^g. At the branch cut < (i^^p^)^ the imaginary 
parts cross and the real parts avoid crossing 

ImA = Ws + ^^^Ki , ^ =5k + 0{k^), 



ReA = - 



(5trDs 
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At f2 = fl^p ^ the crossings of both real and imaginary parts are degenerate 

ReA = -^^^±\J-6k^+0{k'/^), ImA = u;,±^J-Sk^+^-^k+0{k'/^). 
4 4Y Us 4Y LJs 4uJs 

(5.9) 

The evolving eigenvalue branches reconstruct the eigenvalue surfaces shown in 
Fig. 3. In the one-parameter slices of the surfaces the transformation of the eigen- 
value branches from the crossing to the avoided crossing due to variation of the 
parameters O and k occurs after the passage through the exceptional points, where 

the branches touch each other and the eigenvalue surfaces have Whitney's umbrella 
singularities. The surface of the imaginary parts shown in Fig. 3(a) is formed by 
the two Whitney's umbrellas with the handles (branch cuts) glued when they are 
oriented toward each other. This singular surface is known in the physical litera- 
ture on wave propagation in anisotropic media as the double coffee filter (Keck et 
al. 2003, Berry and Dennis 2003). The viaduct singular surface of the real parts 
results from the gluing of the roofs of two Whitney's umbrellas when their han- 
dles are oriented outwards, Fig. 3(b). The double coffee filter singularity is a result 
of the deformation of the MacKay's eigenvalue cone (shown by the dashed lines in 
Fig. 3(a)) by dissipative and non-conservative positional perturbations. The pertur- 
bations foliate the plane ReA = into the viaduct which has self-intersections along 
two branch cuts and an ellipse-shaped arch between the two exceptional points. 
Fig. 3(b). Both types of singular surfaces appear when non-Hcrmitian perturbation 
of Hcrmitian matrices is considered (Kahan 1975, Kirillov et al. 2005). 

Therefore, in a weakly non-Hamiltonian system (2.1) the fundamental quali- 
tative effect of the splitting of the doublets with the definite symplectic (Krein) 
signature is the origination of the double coffee filter of the imaginary parts and 
the viaduct of the real parts. Structural modification of the matrices of dissipative 
and non- conservative positional forces generically does not change the type of the 
surfaces, preserving the exceptional points and the branch cuts. 

6. Unfolding MacKay's cones with mixed signature 

The definite symplectic signature (a/3 > 0) implies £> > and > and thus 
uniquely determines the type of the singular surface for the real and imaginary 
parts of the perturbed eigenvalues. The case of the mixed symplectic signature 
(a/3 < 0) possesses several scenarios for the unfolding of the MacKay's cones by 
the non-Hamiltonian perturbation, because D and N can have different signs. 

When Z) > and > 0, the imaginary parts of the eigenvalues form the double 
coffee filter singular surface whereas the real parts originate the viaduct, Fig. 3. For 
negative D and negative N the type of the surfaces is interchanged: the imaginary 
parts form the viaduct and the real parts originate the double coffee filter, Fig. 4(a). 

Exceptional points are not created for negative values of N/D. In this case the 
eigenvalue surfaces either intersect each other along the branch cut, which projects 
into the line Imc = in the (fi, K)-plane, or do not cross at all. When A^ > the 
surfaces of the imaginary parts ImA(r2, k) cross and the surface ReA(n, k) avoid 
crossing. Fig. 4(b). For N < the surfaces of the imaginary parts are separated 
and that of the real parts cross. Fig. 4(c). 
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Figure 4. Mixed symplectic signature (q/3 < 0): (a) The viaduct ImA(Sl, k) and the double 
coffee filter ReA(r2, k) for D < Q and A'' < 0; (b) the surfaces ImA(r2, k) crossed along 
the branch cut (bold line) and the separated surfaces ReA(r2, k) for D < 0, > 0; (c) 
separated surfaces of imaginary parts and crossed surfaces of real parts for D > 0, A'' < 0. 

7. Example 1. A rotating shaft 

The simplest mechanical systems described by equations (2.1) and (2.2) are some 
two-degrees-of-freedom models of rotating shafts (Kimball 1924, Shieh and Masur 
1968, Nagata and Namachchivaya 1998, Genta 2007). In (Shieh and Masm 1968) 
the shaft is modeled as the mass m which is attached by two springs with the 
stiffness coefficients fci and k2 = ki + k and two dampers with the coefficients /xi 
and /i2 to a coordinate system rotating at constant angular velocity fi. Fig. 5(a). A 
non-conservative positional force I3r acts on the mass. With u and v representing 
the displacements in the direction of the two rotating coordinate axes, respectively, 
the system is governed by the equations (Shieh and Masur 1968) 



mu 
mv 



■ - 



2mQv 
2mflu 



(fc2 - mn'^)v - I3u = 0. 



(7.1) 



In Fig. 5(b) we show a numerically found surface of frequencies for the shaft 
with m = 1 and fci = 4 in the absence of damping and non-conservative forces. The 
surface has four conical singularities corresponding to the doublets ±2z at = 
and to the double zero eigenvalues at the critical speeds = ±2. The cones in 
the subcritical speed range are near-vertically oriented while those at the critical 
speeds are near-horizontal. Consequently, for small stiffness detuning k the system 
is stable in the subcritical speed range and unstable by divergence in the vicinity of 
the critical speeds, where the bubbles of instability in the decay rate plots originate. 

Addition of the non-conservative forces with (3 = 0.2 and damping with /xi =0.1 
and /X2 = 0.2 yields deformation of the conical surfaces with the apexes at r2 = into 
the double coffee-filters. The real parts form the viaduct singular surfaces shown 
in Fig. 5(c) and (d). In the absence of damping (/ii,2 = 0) the gyroscopic system 
with the potential and non-conservative positional forces cannot be asymptotically 
stable in accordance with the theorem of Lakhadanov (1975). It is unstable almost 
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Figure 5. (a) A model of the rotating shaft; (b) four MacKay's cones due to stiffness 
modification {fii — 0, fj,2 = 0, /3 = 0); (c) the viaduct singular surface created by the 
circulatory force only {(3 = 0.2) and (d) by the damping only (/ii — 0.1, fj,2 ~ 0.2). 

everywhere in the space of parameters and can be only marginally stable on the 
set of measure zero in it. This is seen in Fig. 5(c), which shows that the shaft is 
marginally stable at the points of the branch cuts, which form the set of measure 
zero, and unstable at all other points of the parameter plane. 



8. Example 2. A rotating circular string 

Consider a circular string of displacement W{ip,T), radius r, and mass per unit 
length p that rotates with the speed 7 and passes at ip = through a massless 
eyelet generating a constant frictional follower force F on the string (Yang and 
Hutton 1995). The circumferential tension P in the string is constant; the stiffness 
of the spring supporting the eyelet is K and the damping coefficient of the viscous 
damper is D; the velocity of the string in the ip direction has constant value "fr. 
With the non-dimensional variables and parameters 
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Figure 6. (a) The Campbell diagram of the unperturbed rotating string with red dots 
marking the nodes with n = 1; (b) a butterfly distribution (8.8) of the exceptional points 
(open circles) in the subcritical speed range in the (f2, fc)-plane when ^ = G and d — 0.3 
(red open circles correspond to n = 1); (c) projections of the branch cuts (8.9) of the 
coffee filters ImA(ri, k) and the exceptional points for n = 1; (d) projections of the branch 
cuts (8.9) of the viaducts ReA(Sl, k) and the exceptional points for n — 1. 

the substitution of w{ip, t) = u{ip) exp(Ai) into the governing equation and boundary 
conditions yields the boundary eigenvalue problem (Yang and Hutton 1995) 

Lu = \'^u + 2n\u' - {l-n'^)u" = 0, (8.2) 

u{0)-u{2n)^0, u'{0)~u'{27t) = ^-^u{0) + ^^u'{0), (8.3) 

where ' = d^. The boundary eigenvalue problem (8.2) and (8.3) depends on the 
speed of rotation (fi), and damping (d), stiffness (fc), and friction (/i) coefficients. 

We note that the artificialness of the term, corresponding to the non-conservative 
positional forces, in the second of the boundary conditions (8.3), was discussed in 
the literature, see, e.g. (Yang and Hutton 1995, Kirillov 2008). We keep it, however, 
to show how the degeneracy of this operator is seen in the eigenvalue surfaces. 

For d = 0, fc = 0, and fi — the eigenvalue problem (8.2), (8.3) has the 
eigenvalues = m(l +e51), Af„ — im{l + SD.) , where e,6 = ±1 and n,m & Z— {0}. 
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In the (n, IinA)-plane the branches intersect each other at the node {^lo,Uo) with 

^ n — m nm(S — e) ,„ 
^0 = —. , '^0 = ^ (8.4 

mo — ne mo — ne 

where the double eigenvalue Ao = iujq has two linearly independent cigenfunctions 

u1^ = cos{rnp) — £is\a(rup), = cos{mip) — 6i sm(mip) . (8.5) 

Intersections of the branch with n = 1 and e = 1 and the branches with m > and 
(5 < in the subcritical range < 1) are marked in Fig. 6(a) by red dots. 

Taking into account that 5 = — e at all the crossings, excluding (Oq = ±1,wo = 
0) where S = e, we find approximation to the real and imaginary parts of the 
perturbed non-zero double eigenvalues (Kirillov 2008) 



n + m /|c| - Rec 

ReA = —a l^'o ± 

Snnm 



T , 'n — m.^ n + m, , /|c|-|-Rec ,„ 

ImA = wo + e Afi + k±\P-^ , (8.6) 

2 STrnm V 2 ' ^ ^ 

where Afl = fl — fig, and for the complex coefficient c we have 

, 2dLUi) — : 11 ill — m) ( n + m ^ ^ m, — n ,\ f e jn — n 

Imc = k \^ , '--2[e——An + - — /i - d- loq 

167r^nm \ 2 STrnm J \^t^ STrnm 

^ fen — 5m. ^ m-n\'^ fc^ [d{rn + n)ujQf . 

Rec = Af]+- k] +—- \. 2 2 2 ■ 8.7 

Setting Rec = and Imc = we find the coordinates of the projections of the 
exceptional points of the surfaces ReA(f2, k) and ImA(r^, k) onto the (fi, fc)-plane 

ri _L ^ (m + n)(i^a;Q duJo{2efj.nm — d{m — n)u!o) 

ilEP — "0 i 3 , Kep — i — . 

STrnm y'nm{ij?nm + d'^ujl) 2y'nm{^i^nm + d^ujl) 

(8.8) 

As in formulas (5.1) the existence of the exceptional points (8.8) depends on the 
symplectic (Krein) signature of the intersecting branches, i. e. on the sign of nm, 
where n,™ G Z — {0}. In the case of the rotating string all the crossings in the 
subcritical speed range (|n| < 1) have definite Krein signature (nm > 0). For 
those in the supercritical speed range (|r2| > 1) it is mixed with nm < 0. In the 
(f2, K;)-plane the exceptional points are situated on the line Imc = 

/ .2enmu — dujQ(m — n)^^ 

A: = 27r£ n + m , „ '-An. 8.9 

da;o(m2+n2) ^ ^ 

In Fig. 6(b) we show the exceptional points (8.8) of the string passing through 
the eyelet with the damping coefficient d = 0.3. The red open circles correspond 
to the exceptional points born after the splitting of the diabolical crossings with 
n = 1 and e = 1, which arc shown in Fig. 6(a) by the red dots. The exceptional 
points in the {fl, K)-plane are distributed over a butterfly-shaped area, which pre- 
serves its form independently on the number of points involved. In comparison 
with the numerical methods of (Jones 1988, Skorokhodov 2007) our perturbation 
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Figure 7. For d — 0.3 and fi — 0: (a) the double coffee filter singular surface ImA(n, k) in 
the vicinity of the crossing (n = 1, m = 2); (b) the viaduct surface ReA(f2, k) corresponding 
to the crossing {n = l,m — 2); (c) intersecting surfaces lmA(r2, k) in the vicinity of the 
crossing (n = l,m = —2) and (d) the corresponding non-intersecting surfaces KeX{Q,k). 



approach gives efficient explicit and interpretable expressions for the distribution 
of the exceptional points, for the branch cuts, and for the very eigenvalue surfaces. 

In Fig. 6(c) we plot the exceptional points originated after the splitting of the 
diabolical points with n — 1 and e — I together with the projections of the branch 
cuts (8.9) of the double coffee filters lmA(17, k), which are shown by the bold lines. 
The corresponding projections of the branch cuts (8.9) of the viaducts ReA(il, k) 
are presented in Fig. 6(d). Only exceptional points originated after the perturbation 
of the doublets with fio = are situated on the f2-axis. This explains why damping 
creates a perfect bubble of instability for the doublets with m = n and imperfect 
ones for the diabolical points with m ^ n (Yang and Hutton 1995, Kirillov 2008). 
Approximations (8.6) to the eigenvalue surfaces of a string with /i = and d — 0.3 
are presented in Fig. 7 for different values of n, m, e, and 5. The smaller inclusions in 
Fig. 7 show the cross-sections of the surfaces by the plane fc = for the convenience 
of comparing with the numerical data of Yang and Hutton (1995). The results 
shown in Fig. 7 are in qualitative agreement with the developed theory for the 
equations (2.1) and (2.2) and perfectly agree with the numerical modeling. 

In Fig. 8 we show in the complex plane the parent diabolical points (red dots) 
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Figure 8. Exceptional points (open circles) with the parent diabolical points with n — 1, 
e = 1 in the complex plane when d — 0.3 and the trajectories \{Q,) for k = 0.05. 

and the corresponding exceptional points (open circles) whose locations are 



47r n + m 



2nm ^ d n — m 



(8.10) 



In the engineering literature it was observed that the exceptional points (strong 
modal resonances (Dobson et al. 2001)) are precursors to flutter instability because 
of their strong influence on the movement of eigenvalues in the complex plane. Fig. 8 
demonstrates the approximation of the 'dynamics' of eigenvalues in the vicinity 
of the exceptional points, calculated by the formulas (8.6), which is in a good 
qualitative agreement with the known numerical results (Sinou et al. 2006). 

Finally, we notice that pure non-conservative positional perturbation [d — 0) 
causes degeneration of the eigenvalue surfaces. Indeed, the line (8.9) reduces to 
= f2o and the two exceptional points merge into one at ^ep = and kep — 0. 
As a consequence, the central arch of the viaduct and the branch cut of the double 
coffee filter shrink to a single point. At this exceptional point the angle of crossing 
of the surfaces is zero in agreement with (Yang and Hutton 1995, Kirillov 2008). 
This degeneration visualizes the artificialness of the term related to the friction 
force in (8.7) that was already pointed out in the literature by physical arguments. 



9. Conclusion 

We found that in a weakly anisotropic rotor system (2.1) the branches of the Camp- 
bell diagram and the decay rate plots in the subcritical speed range are the cross- 
sections of the two companion singular eigenvalue surfaces. The double coffee filter 
and the viaduct are the imaginary and the real part of the unfolding of any double 
pure imaginary semi-simple eigenvalue at the crossing of the Campbell diagram with 
the definite symplectic (Krein) signature. Generically, the structure of the perturb- 
ing matrices determines only the details of the geometry of the surfaces, such as the 
coordinates of the exceptional points and the spacial orientation of the branch cuts. 
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It does not yield the qualitative changes irrespective of whether dissipative and cir- 
culatory perturbations are applied separately or in a mixture. The two eigenvalue 
surfaces found unite seeming different problems on friction-induced instabilities in 
rotating elastic continua, because their existence does not depend on the specific 
model of the rotor-stator interaction and is dictated by the symplectic signature 
of the eigenvalues of the isotropic rotor and by the non-conservative nature of the 
forces originated at the frictional contact. The double coffee filter singularity and its 
viaduct companion are true symbols of instabilities causing the wine glass to sing 
and the brake to squeal that connect these phenomena of the wave propagation in 
rotating continua with the physics of non-Hermitian singularities associated with 
the wave propagation in stationary anisotropic chiral media (Berry 2004). 
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